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ABSTRACT 



We present VLT/UVES spectroscopy of the quasar Q0841+129, whose spectrum shows a proximate damped Ly-a 
(PDLA) absorber at 2=2.47621 and a proximate sub-DLA at 2=2.50620, both lying close in redshift to the QSO itself 
at 2:cm=2.49510±0. 00003. This fortuitous arrangement, with the sub-DLA acting as a filter that hardens the QSO's 
ionizing radiation field, allows us to model the ionization level in the foreground PDLA, and provides an interesting 
case-study on the origin of the high-ion absorption lines Si iv, C iv, and O vi in DLAs. The high ions in the PDLA show 
at least five components spanning a total velocity extent of ~160kms~^, whereas the low ions exist predominantly 
in a single component spanning just 30kms~^. We examine various models for the origin of the high ions. Both 
photoionization and turbulent mixing layer models are fairly successful at reproducing the observed ionic ratios after 
correcting for the non-solar relative abundance pattern, though neither model can explain all five components. We show 
that the turbulent mixing layer model, in which the high ions trace the interfaces between the cool PDLA gas and a 
hotter phase of shock-heated plasma, can explain the average high-ion ratios measured in a larger sample of 12 DLAs. 
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1. Introduction 

In the process of galaxy assembly and evolution, interstel- 
lar and circumgalactic gas clouds play a key role, represent- 
ing both the fuel for future star formation and the metal- 
enriched by-products of past star formation. Absorption- 
line spectroscopy is a powerful way to study these gas 
reservoirs, to characterize their kinematic and chemical 
properties, and to investigate their link to star formation. 
Of particular interest for probing interstellar gas are the 
damped Ly-a systems (DLAs), the highest column density 
absorbers defined by log7V(Hi)>20.3 (Wolfe et al. 198(|. 
Their high neutral gas content ensures they select over- 
dense regions, including galaxies and their halos. DLAs con- 
tain the majo rity of the neutral gas in the Universe in the 
range z=0-5 (jWolfe et al.l 120051 ) . and serve as luminosity- 
independent indicators of cosmic chemical enrichment. 

Despite being known for their neutrality, DLAs fre- 
quently show high-ion absorption overlapping in velocity 
space with the neutral lines yet tracing a separate phase 
(or phases) of gas. The first high ions detected in DLAs 
were Civ and Si iv ( IWolfe et al.l ll994t iLu et al.l[T995l Il996t 
I Wolfe fc Prochaskall2000ai rbl). which almost always show a 
different (though overlapping) component structure than 
the low ions. We have recently been surveying high-ion ab- 



* Based on observations taken under proposal IDs 65.0- 
0063(B), 70.B-0258(A), and 383.A-0376(A) with the Uhraviolet 
and Visual Echelle Spectrograph (UVES) on the Very Large 
Telescope (VLT) Unit 2 (Kueyen) operated by the European 
Southern Observatory (ESO) at Paranal, Chile. 



sorption in DLAs at z=2-4 at high resolution using the 
VLT/UVES spectrograph, finding the high ions are ubiq- 
uitous. In data with S/N ^30 per pixel, we report detection 
rates of > 34% for Ovi (iFox et al.l l2007al). 100 % for CiV 
and Si IV (jFox et al.ll2007bD . and 13% for Nv (jFox et al.l 
,2_009) . High-ion absorption is also see n in a composite SDS S 
sp ectrum of several h undred DLAs (jRahmani et al.ir2010[ ). 
In iFox et all (|2007aD we favored the interpretation that 
the O VI absorption components in DLAs trace warm-hot 
(T ^ 10^ K) coUisionally-ionized plasma, rather than cool 
(T ~ 10'* K) photoionized plasma, based on two obser- 
vational findings: first, a tendency for the individual O vi 
components in DLAs to be broader than the CiV compo- 
nents (suggesting higher temperatures in the O vi phase) , 
and second, the failure of photoionization models to repro- 
duce the Ovi column densities measured in a particular 
DLA studied in detail (the z^bs=2.076 DLA toward Q2206- 
199) without recour se to unreasonably large path lengths. 
iLehner et al] (|2008[ ) investigated the Zabs=2.377 DLA to- 
ward J1211-I-0422, and similarly found that photoionization 
cannot reproduce the observed O Vl column, although they 
also found that single-phase collisional ionization models 
(equilibrium or non-equilibrium) cannot explain the high- 
ion column densities in that system. 

The origin of Ovi absorption in DLAs demands fur- 
ther investigation, since the observational characterization 
of warm- hot gas in galaxy halos at redshifts 2-3 would pro- 
vide useful constraints on models of galaxy formation and 
evolution. Unfortunately, observations of O vi absorption at 
these redshifts are not straightforward, since the O vi reso- 
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nance doublet (AA1031,1037) lies in the Ly-a forest and so 
often suffers from severe contamination from foreground H i 
absorbers (e.g. Buries & Tvtler 1996; Frank et al. 2010). 
With a large enough DLA sample, one occasionally comes 
across a case where the Ovi doublet lines are only par- 
tially (rather than completely) blended, and where informa- 
tion on portions of the O VI l ine profiles c an be recovere d 
(JFox et al. 2007 J: iLehner et al. 2008: Simcoe et al.ll2002D . 
Such a case (the DLA at z^bs='2A7621 toward Q0841+129) 
is discussed in this paper. This system has another advan- 
tage: although it lies at only 1630 km s^^ from the back- 
ground QSO (classifying it as a proximate DLA or PDLA) , 
it is shielded from the QSO's ionizing radiation by a second 
absorber, a sub-DLA at Zabs=2. 50620, very close in redshift 
to the QSO at Zom=2.49510±0.00003. The H i column in the 
sub-DLA, log7V(Hi)=19.00±0.15, is high enough that the 
system is optically thick at the Lyman Limit, but not at 
high photon energies {E^75eV, i.e. A ^165 A). Therefore 
the sub-DLA acts to harden the ionizing radiation field in- 
cident on the PDLA. 

PDLAs are defined as those DLAs lying within 
5 000 km s~^ (or sometim es 3 OOP km s~^) of the r e dshift 
of the background QSO (lEllison et all l2002l l20Tol l20lH 
iRussell et a"Lll2006l : lRbret al.ll2007t iProchaska et al.ll2008D . 



They are frequently excluded from studies of intervening 
DLAs, due to their potentially unusual ionization condi- 
tions and environmental effects (such as galaxy clustering 
near quasars). However, the fact that the (enhanced) ioniz- 
ing radiation field in PDLAs is calculable can be turned to 
our advantage, and allows them to be used as case-studies 
for the origin of high-ionization plasma in DLAs. That is 
the approach followed in this paper. 

2. Observations and Data Handling 

Q0841-H29 has been observed with VLT/UVES 
(jDekker et al.l |2000[ ) under three ESO programs, each 
time with different dichroic settings and wavelength 
coverage. These data were taken in 2000 (Dic2 412-1-860 
setting, program P65.O-0063), in 2002 (Dic2 390-f860, 
70.B-0258), and in 2009 (Did 390-^580, 383.A-0376), ah 
with a 1.0" slit (except the red-arm P70 observations, 
which used a 0.9" slit). We downloaded these data and 
their associated calibration files from the E SO archive, and 
redu ced them using the UVES pipeline (Ba llester et al.l 
I200Q!) in MIDAS, co-adding the individual exposures after 
weighting by signal-to-noise. This produced a final spec- 
trum with 6.6kms~^ resolution (FWHM) in «2kms~^ 
pixels covering the range 3 286 to 10 429 A (with a gap 
between 5 756 and 5 837 A). 

Three DLAs are present along the Q08 41-I-129 sight 
line, at Zabs=l-8639, 2.3745, and 2.47621 (jLedoux et al.l 
|2006| ) in the vacuum heliocentric frame. In addition, a 
sub-DLA is detected at Zabs=2. 50620, with all the hall- 
marks of an intrinsic absorber (see §6). The proxim- 
ity (1630 kms-i) of the z=2.47621 DLA to the QSO 
redshift Zem=2.49510±0.00003 (derived in §4) classifies 
it as a PDLA. The system has lo giV(Hi)=20.80±0.1 
and [Zn/H] = [Znii/Hi]=-1.60±0.10 ( Ledoux et al.l 120061) . 
a good indication of the true metallic ity since Zn is typ i- 
cally undepleted onto dust grains (e.g. lPettini et al.lll997D . 
In a WMAP 7-yea r cosmology (iifo= 71 kms""'^ Mpc""'^, 
1^M=0.27, n^=0.73■. [Jarosik et al.ll201lh . the Hubble pa- 
rameter at z=2.5 is iJ(2.5)=250 kms"-'^ Mpc"-'^ (using the 



relation given in iPhillipps et al.l l2002f) . Therefore, assum- 
ing the relative motion between the QSO and PDLA is 
purely due to Hubble flow (i.e., assuming peculiar veloci- 
ties are negligible, though see Ellison et al. 2010), we find 
the PDLA lies at «6.5 proper Mpc from the quasar. 

3. H I column densities in proximate DLA and 
sub-DLA 

In Figure 1 we show the Ly-a profiles of the two prox- 
imate absorbers. Ly-a absorption from the sub-DLA at 
-2abs=2. 50620 hes in the red wing of the Ly-a absorption 
from the DLA at Zabs=2. 47621. The sub-DLA is far-enough 
separated from the DLA that its own damping wings are 
well-defined and can be fit to determine its Hi column 
density. These wings are visible in Figure 1 at observed 
wavelengths of 4250-4260 and 4265-4270A. Voigt-profile 
fits to the two absorbers yield log Ar(H i)==20.80±0.10 in the 
PDLA and logiV(Hi)=19.00±0.15 in the proximate sub- 
DLA. 




4180 
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Fig. 1. Ly-a profile in the two proximate absorbers toward 
Q0841-hl29. The red line shows the best-fit Voigt profile 
with log7V(Hi)=20.80±0.10 in the DLA at Zabs=2.47621 
(centre of the figure) and logiV(Hi)=19.00±0.15 in the 
sub-DLA at 2:abs=2. 50620 (toward the right), la errors are 
shown with blue dashed lines. The red dotted line shows 
the fit to the DLA only 



4. Emission-line redshift and specific luminosity of 
Q0841+129 

No precise measurement of the redshift of Q0841-I-129 cur- 
rently exists in the literature. We used a flux-bisecting al- 
gorithm to measure the observed emission-line wavelengths 
of Ly-/3 1025.722, Oi 1304.460 (multiplet), Siiv 1393.755, 
Si IV 1402.770, Nv 1240.14 (doublet). Civ 1548.195, Heii 
1640.42, Cm] 1908.73, and Mgii 2798.74. A weighted av- 
erage of the redshifts determined from each of these nine 
lines gives Zom=2.49510±0. 00003, the value we adopt for 
the rest of this paper. Although it is common practice to 
apply velocity offsets to each line before averaging (e.g. 
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IVanden Berk et al]|2001h . we found no need to do this. The 
Mr II red shift (considered a good measure of the true value; 
iHewett fc Wild 2010) is 2.49512±0.00001. 

We derived the specific lu minosity of Q0841-I-129 at the 
Lyman Limit, Lqi2, following iGuimaraes et all ()2007l ). The 
method uses the QSO's observed ^-magnitude (18.5) and 
emission-line redshift Zem, assumes an intrinsic QSO spec- 
tral slope of —0.5 (where F^oci^"), and corrects for Galactic 
extinction using the value E( B — y) =0.048 found t oward 
Q0841-I-129 from the maps of ISchlegel elall (|1998[ ). This 
gives i9i2=8.86xlO^°ergs~^Hz~^. Assuming the QSO- 
PDLA distance is Hubble-Flow dominated (6.5 Mpc), we 
derive a specific flux of QSO radiation at 912 A at the 
PDLA of F^l° ', 



Table 1. Apparent optical depth measurements of low- 
ion absorption in the Q0841-I-129 PDLA 



1.7 X lO^^ergcm ^s ^ Hz 



a factor of 



«6 lower than the calculated fiux of the UV background 
(UVB) at this redshift fjjy^ » lO-'^^ergcm"^ s ~^ Hz ~ ^ 
JHaardt fc Madaul 120011: iFaucher-Giguere et all l2009f PI. 
Based on this calculation of Fgi2, the PDLA is outside the 
proximity zone of the quasar, although for higher photon 
energies, we find the PDLA does lie inside the proximity 
zone (see §6). Despite these complications in how to de- 
fine proximity, we still refer to the absorber as a PDLA, 
following convention. 

5. Absorption in PDLA 

In Figure 2 we show the normalized profiles of 36 ab- 
sorption lines of both high and low ions in the PDLA. 
The low-ion profiles show a simple component structure; 
Sii A1259 shows a single unsaturated component spanning 
Ad9o=30 kms~^ (where Awgo is the width containi ng th e 
central 90% of the optical depth; Prochaska &: WolfQil997l) . 
An additional weak low-ion component is seen in G ii and 
Sill near 50kms~^ in the PDLA rest-frame. Integrated col- 
umn densities for the low ions were me asured using the ap- 
paren t optical depth (ADD) technique (Sav age fc SembachI 
Il991h . using atomic data from sMorton (20031) . and are pre- 
sented in Table 1. 

In the high ions, multi-component absorption is ob- 
served in the PDLA in Siiv, Giv, and Ovi, though the 
Ovi profiles are partly blended. Nv is not detected. Voigt- 
profile fits were conducted to the high-ion absorption using 
the VPFIT software package0. These fits are over-plotted 
on Figure 2 as red lines and detailed in Table 2, where we 
also present AOD measurements of the high ions. We chose 
not to tie the GiV and SilV component centroids during 
the fit process, since several differences between the pro- 
files suggest the two ions are not fully co-spatial. 

The high-ion component structure in the PDLA differs 
from the low-ion structure. This is illustrated in Figure 
3, where the apparent column density profiles of several 
high ions and Fell A1608 (chosen as an unsaturated low- 
ion line) are compared. The high ions are contained in 
several components extending over wl60kms~^, in sharp 
contrast to the narrowness of the Fe ii line covering just 
«30kms~^. Based on the GiV profile we identify five high- 
ion components in the PDLA, named as follows: Gl cen- 



A recent measurement of the UVB at 2=2.73 using 
the proximity effect i s F^^^ ^ 4 x 10"^^ ergcm"^ s"^ Hz"^ 
IjDall'Aglio et al.] 120081 1 . several times lower than the Haardt & 
Madau (2001) value, but still higher than the flux incident on 
the PDLA from the background QSO. 

2 



Line 



(kms-^) 



(mA) (Afa 



% 



logA^^ 



Ni 1134.2 


-30 


30 


13±3 


13.96±0.03 


Ni 1134.4 


-30 


30 


26±3 


13.99±0.01 


Nil 1083 


-30 


30 


44±3 


13.69±0.01* 


Sii 1250 


-30 


30 


22±3 


14.51±0.02 


Sii 1253 


-30 


30 


49±3 


14.57±0.01 


Sii 1259 


-30 


30 


59±3 


14.53±0.01 


Cii 1036 


-30 


30 


180±3 


>14.72 


C II 1334 


-30 


30 


243±3 


>14.65 


Cm 977 


-30 


30 


200±3 


>14.17 


An 1048 


-30 


30 


24±3 


13.05±0.02* 


Sin 1260 


-30 


30 


241±3 


>13.85 


Sin 1193 


-30 


30 


207±3 


>14.06 


Sin 1190 


-30 


30 


179±4 


>14.27 


Sin 1526 


-30 


30 


220±3 


>14.52 


Sin 1304 


-30 


30 


173±3 


>14.59 


Sin 1808 


-30 


30 


54±3 


14.99±0.01* 


Mg II 2796 


-30 


30 


505±3 


>13.60 


Mgii 2803 


-30 


30 


494±4 


>13.92 


Zn II 2026 


-30 


30 


12±3 


11.83±0.04 


Cr II 2056 


-15 


15 


21±3 


12.76±0.02 


Cr II 2062 


-15 


15 


20±3 


12.87±0.02 


Cr II 2066 


-15 


15 


10±3 


12.74±0.03 


Fell 1144 


-30 


30 


119±3 


>14.34 


Fe II 1608 


-30 


30 


161±3 


>14.51 


Felll 1122 


- 5 


22 


21±3 


13.58±0.01* 


Aim 1854 


-20 


20 


56±3 


12.59±0.01 


Aim 1862 


-20 


20 


30±3 


12.60±0.01 


Sim 1206 


-30 


30 


249±3 


>13.75 



See |http://www.ast. cam. ac.uk/^rfc/vpfit. html 



" Velocity range used in AOD integration, relative to 
2abs=2.47621. The default is -30 to 30kms"\ but this is 
reduced when necessary to avoid blends. 

^ Rest-frame equivalent width, with error that includes 
statistical and continuum-placement contributions combined 
in quadrature with a 2.5 mA error accounting for the uncer- 
tainties in the velocity integration range. 

'^ Apparent column density. Lower limit given if line is sat- 
urated. Errors include statistical and continuum-placement 
contributions combined in quadrature with a 0.01 dex error 
accounting for the uncertainties in the velocity integration 
range. Entries marked * are potentially blended and should 
be treated with caution. 



tered at — 60kms^^, G2 at — 15kms^^, G3 at 15kms^^, 
G4 at 55kms~^, and G5 at 102kms~^ (all defined rela- 
tive to Zabs)- Although GiV and Si IV show similarities in 
their component structure, there are clear differences: Si IV 
is detected in G2, G3, and G4, but not in Gl and G5. The 
relative strength of Si iv to G iv changes substantially be- 
tween G2, G3, and G4. Furthermore, G2 and G4 are broader 
in Giv than in Siiv (see &- values reported in Table 2). 

Partial blending complicates the measurement of O vi. 
We report an Ovi detection in the range —100 to Okms^"'^ 
(corresponding to components Gl and G2; see Figure 3), 
because the two lines of the O vi doublet show largely con- 
sistent optical depth profiles at these velocities (although 
small differences are observed at —70 to — 40kms^^). 
Furthermore, the integrated AOD column densities of the 
two O VI lines made over the interval — 100 to km s~ ^ agree 
to within 0.04 dex (see Table 2). Apparent features in Ovi 
A1031 and O vi A1037 follow G iv components G3 and G4, 
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PDLA at z = 2. 47621 
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log N=13.04 
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log N=14.81 
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log N=11.73 
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log N=12.82 
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log N=12.82 
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log N=14.65l 
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log N=12.59 
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log N=12.59 
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1 K~-^f'A-^/^AUU/]jy'fij^-\yuV''-h 



Si IV 1402 

■nwrrt 



log N = 13.63 



Si III 1206 




log N<12.57 





1 \r^v^J'yY\vMiMt'i"^'^j--^yi-yi^f'J 



N V 



1242 
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Fig. 2. Normalized VLT/UVES absorption-line profiles in the PDLA at z=2.47621 toward Q0841+129. Red lines show 
Voigt-profile fits, with tick marks showing component centers. Fits were not attempted for Sim, Cm, or Mgii due to 
extreme saturation. Grey shaded regions indicate blends. The total column density obtained from component fitting 
(summed over components) is annotated on each panel. 



respectively, but in each case the data are noisy, and the other O vi line is blended, so we cannot confirm the detec- 
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Table 2. Voigt-profile fits and AOD measurements of high-ion absorption in the PDLA toward Q0841+129 



Ion Comp. vq^ b fog TV iogiVtoi^ 
(knis~^) (kms~^) (JV in cm~^) (JV in cm" 



iog7Va(strong;[^ 
(JV in cm~^) 



log A''a (weak)'' 
(JVin I 



') ( 



V- ,v+_ 
kms^^l 



Si IV 


CI 






<11.72 


13.63±0.01 


13.60±0.01 


13.63±0.01 


-50,80 




C2 


-13±1 


16.7±0.4 


13.34±0.01 




13.28±0.01 


13.28±0.01 


-100, 




C3 


15±1 


9.9±0.6 


12.88±0.02 












C4 


51±1 


12.5±0.3 


13.06±0.01 












C5 






<11.86 










Civ 


CI 


-57±3 


19.2±3.2 


12.78±0.13 


13.89±0.03 


13.89±0.01 


(blended) 


-100,120 




C2 


-12±1 


27.3±3.1 


13.53±0.04 




13.49±0.01 


13.49±0.01 


-100, 




C3 


17±1 


8.8±1.8 


12.69±0.14 












C4 


55±1 


15.8±0.3 


13.48±0.01 












C5 


102±1 


9.4±1.5 


12.40±0.05 










Ovi^ 


CI 
C2 


-55±3 

- 8±5 


25.0±2.5 
24.3±17.5 


13.87±0.06 
13.39±0.32 


13.99±0.11 


13.97±0.01 


13.93±0.01 


-100, 


Nv^ 


CI 










<12.57 


<12.86 


-100, 





















" Velocity centroid relative to 2abs=2. 47621. 
Total VPFIT column density, summed over all components, with errors added in quadrature. 

'^ Apparent column density measured from strong and weak doublet lines in velocity interval indicated. 

** Velocity range used in AOD integration, relative to Zabs. 

"^ Voigt fit for Ovi only possible in components CI and C2. Ovi AOD measurements made over velocity region where the 
profiles of the two doublet lines overlap. The Civ, Siiv, and Nv measurements are repeated over this interval so the ionic ratios 
can be derived. 

^ No detection in either Nv line. Upper limits are 3a. 



VI I03I 
VI 1037 
Si IV 1393 (x3) 



• C IV 154H (xa) 
O C IV 1550 (x2) 

• Fe II 1608 (x0.15) 

I 




Q0841+129 
2=2.47621 

.<3 







V (km s ') 



Fig. 3. Top panel: Apparent column density profiles of high-ion and Fell absorption in the Zabs=2. 47621 PDLA. Several 
of the profiles have been scaled by the factors shown in the legend, to allow inter-comparison with the other lines. The 
O VI profile is only shown at velocities where the two O vi Unes overlap; upper limits for O vi are shown at other velocities, 
using whichever of the two doublet lines gives the stronger constraint. Middle panel: C iv/0 vi apparent-column-density 
ratio as a function of velocity; in pixels with upper limits on O VI, lower limits on the ratio are given. Lower panel: 
Civ/Siiv apparent-column-density ratio as a function of velocity. The mean (solid lines) and standard deviation (dashed 
lines) of the logarithmic ratio in each component are shown in red. 



tion in those components. We therefore only report upper 
limits on Ovi in C3, C4, and C5. For each of these three 



components the upper limit was derived from whichever of 
the two O VI lines gave the stronger constraint. These up- 
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per limits are used to constrain the high-ion ratios in each 
component. 

The probabihty that the O vi detections we report in 
CI and C2 are genuine can be roughly estimated using the 
known line density of the L y-g forest, a majo r contami- 
nant. The data presented by iKim et al.l ()1997[ ) indicate a 
mean separation of w500kms~^ for Ly-a forest lines with 
7V(Hi)=13. 1-14.3 at z=2.5, indicating that the probabil- 
ity of a forest line falling in a 20kms~^ interval around a 
given redshift is Pbiond ~ 20/500 w 0.04 (ignoring cluster- 
ing) . The probability of two forest lines falling at that red- 
shift, one in each line of Ovi, is ^biond ~ 0.0016. However, 
the observed density of blends is far higher than one per 
500kms~^, indicating that other contaminants (metal lines 
and higher-order Lyman lines) are also present. We mini- 
mize the chance of these contaminants mimicking O VI by 
relying on intervals where both O vi lines agree. 



6. Absorption in proximate sub-DLA 

Absorption-line profiles from the proximate sub-DLA at 
-2abs=2. 50620 are shown in Figure 4. This system has 
logiV(Hi)=19.00±0.15 and shows strong high-ion absorp- 
tion and weak low-ion absorption. In particular, it ex- 
hibits strong multi-component N v absorption, which is 
usually absent from intervening DLAs, multi-component 
Svi A944 absorption, also unknown in intervening DLAs, 
and fully saturated O vi absorption. It has an unusually 
low ratio [Oi/IIi]=— 2.43±0.05, indicating unusual ion- 
ization conditions given that intri nsic sub-DLA metal- 
hcities are rarely that low (e.g iKulkarni et al.l l2007t 
iDessauges-Zavadskv et al.ll2009f) . All these properties iden- 
tify the sub-DLA as an intrinsic system arising close to the 
central engine of the AGN. The intrinsic nature of the sub- 
DLA supports our interpretation that it lies between the 
PDLA and the QSO, as is suggested (but not necessarily 
required) by the sub-DLA having a higher redshift than the 
PDLA. 

The intrinsic sub-DLA could be used to study the phys- 
ical conditions and abundances in the gas immediately sur- 
rounding the AGN. However, in our case, the system is 
of principal interest as a blocking filter, serving to attenu- 
ate the background QSO radiation incident upon the fore- 
ground PDLA. The attenuating effect of the sub-DLA on 
the QSO's ionizing radiation can be calculated from the 
following equations, which relate the Hi and Hen optical 
depths to the H i and He ii columns and wavelengths: 
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where crHi=6.3xlO~^® cm^ and crHcii=crHi/4 are the 
photoionization cross sections for Hi and Hen, respec- 
tively. These expressions are evaluated using the mea- 
sured sub-DLA Hi column log A^(Hi)=19.00, assuming 
He/H=0.0823 by number (equivalent to a He mass frac- 
tion of 0.24), and assuming all the He is in the form of 
Hen, giving logiV(Heii)=17.9. The total optical depth is 
then T=THi+THeii (we do not consider metal-line opacity). 
The effect of the H i and He n opacity from the sub-DLA 
is to remove photons with energies 13.6 < E ^ 75 eV from 
the radiation field, with the Hi providing the opacity at 
13.6 < E < 54eV and the Hen at 54 < S < 75 eV. Thus 



Fig. 4. Normalized VLT/UVES absorption-line profiles of 
the proximate sub-DLA toward Q0841-I-129. The velocity 
zero-point shown is defined by the strongest N v absorption 
component (2=2.50552), slightly lower than the redshift de- 
fined by the low ions (2=2.50620). The presence of strong 
Svi and Nv absorption, and fully saturated Ovi absorp- 
tion identifies the sub-DLA as intrinsic to the quasar. 



the sub-DLA hardens the ionizing radiation field incident 
on the PDLA, by preferentially transmitting photons at 
i?<; 75 eV. This is shown in Figure 5, where the total radia- 
tion field Fi, incident on the PDLA is derived. The radiation 
field is calculated as the sum of: (1) the quasar spectrum 
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attenuated by the sub-DLA, assuming an intrinsic power- 
law spectrum F^, oc v" with q;=— 0.5, a normahzation from 
the value of Fqi2 derived in §4, and attenuated using the 
sub-DLA optical depths described above, and (2) the extra- 
galactic U VB at 2:= 2. 5, which is built int o CLOUDY and 
is based on lHaardt fc Madaul ()l996ll200lt) . We do not con- 
sider internal sources of radiation within the PDLA (e.g. 
any star-forming regions), nor do we consider any attenu- 
ation of the radiation field from the PDLA itself. Figure 
5 shows that the UVB is stronger than the QSO radia- 
tion field at energies E ^75 eV, but above this energy, the 
QSO radiation dominates. The QSO radiation is therefore 
important to include when investigating the origin of O VI 
(which requires 114eV for its creation). The combined ra- 
diation field is used as an input for the CLOUDY modelling 
described in §8. 
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Fig. 5. Spectra of the ionizing radiation fields incident upon 
the PDLA: the quasar spectrum attenuated by the sub- 
DLA (green), the extragalactic UV background (UVB) at 
z=2.5 (blue), and the sum of the two (orange). Specific 
flux is plotted against frequency, with the ionization edges 
of various ions annotated on the plot. The derivation of 
the spectra is discussed in §6. Note how the QSO spectrum 
dominates at extreme UV energies {E > 75 eV) . 



7. PDLA ionization level: the low-ion phase 

The ionization level in the low-ion phase of DLAs 
can be diagnosed with the Arl/Sili, A r l/Sli , 
Feiii/Feii, and Nii/N i ratios (Prochaska et al.l 120021: 
iHowk fc SembachI Il999t IvTadilo et al. 2001). In the 
PDLA, we measure log A^(Ari)/A^(Siii)=-1.77±0.15, 
log A^ (Ar i)/N{S ii)=-1.47±0.03, 

logA''(Feiii)/Ar(Feii)=-1.17±0.05 (or -1.02±0.05 
when considering the strongest component only), and 
log7V(Nii)/A^(Ni)=-0.31±0.07, using the integrated 
VPFIT column densities. Each of these four ratios indi- 
cates the low-ion gas is predominantly (but not completely) 



neutral with a hydroge n ionization fraction x(IIi)=0.1-0.5 
(jProchaska et al.l |2002[ ): a value a;(Hi)>0.5 is ruled out 
since that would lead to logA^(Nii)/Af(Ni)>-0.2 and 
logAf(Feiii)/Af(Feii)>-1.0; a value x(Hi)<0.1 is ex- 
cluded since that would lead to log Af(Ari)/Ar(Sili)>-0.2 
(ISofia fc Jenkins! I1998D and log Af(Ari)/A^(Sii)>-0.2 
(neg lecting any poss i ble du st or nucleosynthesis effects). 
The iProchaska et al.l (J2002t ) model predictions were based 
on CLOUDY calculations that used the UV background 
as the incident radiation field. They are applicable to the 
PDLA since, for the energies of interest for the low ions 
{E ^ 50 eV), the radiation field incident on the PDLA 
is UV-background dominated (Figure 5). Note that the 
measured Felll/Feii ratio in the PDLA is technically an 
upper limit as Feiii A1122 is detected at low significance 
close to several blends, and there is an sslO kms~^ offset 
between the Fen and Feiii velocity centroids. However, if 
the feature we fit as Feiii is a blend, then the true Feiii 
column and Feiii/Feii ratio will be lower, thus strength- 
ening the conclusion that the low-ion gas is predominantly 
neutral. Likewise, any contamination of N ii A1083 would 
lower the true Nil column and the Nii/Ni ratio, again 
reinforcing the finding that the low-ion gas has x(Hi)<0.5. 

8. PDLA ionization level: the high-ion phase 

In this section we use the observed high-ion column den- 
sities and their ratios to test various models of the ori- 
gin of the high-ion plasma. The C iv/Si iv, C lv/0 VI, and 
Nv/Ovi ratios measured by Voigt profile fitting in each 
PDLA component are given in Table 3, with limits given 
where appropriate. We also list the average PDLA high- 
ion ratios integrated over velocity, and for comparison, the 
high-ion ratios measured in a range of other galaxy-halo 
environments. Considerable variation in the ratios is seen 
between components, indicating variable ionization condi- 
tions; the Civ/Siiv ratio varies from 0.65±0.25 in C3 to 
>11.5 in CI. The Civ/Ovi ratio varies by over a factor 
of ten between the two components where we detect O vi, 
from 0.08±0.03 in CI to 1.4±1.0 in C2. 



8.1. Photoionization models 

The first high-ion origin scenario we examine is pho- 
toionization (PI), which is often considered the ioniza- 
tion mechanism for the C i v and Si iv absorp tion in DLAs 
(iLu et al.|[l995t iMo fc Mir alda-Escudc 1996|; iMaller et all 
120031: IWolfe et al.ll2005l : lLehner et al.ll200 8). but which has 



been una ble to explain th e D LA O Vi observation s pre- 
sented bv lFox et all ()2007aD and lLehner et al.l (|2008[ ). This 
model is of particular interest in the case of the PDLA 
given the enhanced EUV radiation field incident on the ab- 
sorber. We used th e CLOUDY PI code (last described in 
iFerland et al.|[l998r ) to generate models for PDLA compo- 
nents CI and C2, the only two with Ovi measurements. 
The CLOUDY models assume the gas has uniform density 
and plane-parallel geometry, and use the two-component 
(QSO+UVB) radiation field shown in Figure 5. The model 
uses the measured metallicity [Zn/II]=— 1.60 and a DLA 
relative abundance pattern with [C/0]=— 0.5, [Si/O]=0, 
and [N/O1=-1.0 , typical for DLAs with -2 ^ [0/H] <> 
-1 (iPettini et al. 2008; P etitiean et al.ll2008D. S olar abun- 



dances were taken from lAsplund et al.l ()2009l ). No con- 



straint on A^(Hi) in each high-ion component is avail- 
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Table 3. Component high-ion ratios: PDLA vs averages 
in other galaxy-halo environments 



Component (velocity) 



CI ( -57kms 
C2 ( -12kms" 
C3 ( 17kms"^) 
C4 ( 55kms*^) 
C5 ( 102 km s"^ 
PDLA average,^ 



iV(Clv) a 
N(Si IV)- 



> 11.5 
1.54±0.17 
0.65±0.25 
2.64±0.07 

> 3.4 
1.95±0.09 



NjCn) 
N(0\ 



0.08±0.03 

1.4±1.0 

>0.25 

>0.56 

>0.06 

0.35±0.02 



A/(Nv) 
iV(Ovi) 



<0.049 
<0.194 



<0.042 



DLA average_^ 
Milky Way disk average_^ 
Milky Way halo average_^ 
LMC average;^ 



5.37±2.20 
3.3±2.5 

4.19±1.47 
1.5-2.5 



0.56±0.28 
0.60±0.47 
0.71±0.36 
<0.09-0.59 



" Component ratios derived from VPFIT columns, with 
limits given for components with Siiv and Nv non- 
detections, or potential Ovi blending. 

' PDLA average ratios derived from AOD integrations. 
The C iv/Silv ratio is measured over the full interval —100 to 
120kms~^. The Civ/Ovi and Nv/O vi ratios are measured 
over the (unblended) velocity interval —100 to Okms~^. 

"^ Average and standard deviation of ratio among 12 DLAs 
with Ovi absorption at 2=2-3 (iFox et all l2007al l. The 
N v/O VI average is an upper limit since N v is only detected 
in 3 of these 1 2 DLAs. 

'^ Taken from lLehner et all (|201ll '). Based on measurements 
made along 38 Galactic disk sight lines. 

"^ Taken from IWakker et al.l (|201ll ). Based on measure- 
ments made along 58 extragalactic sight lines. 

^ Range obser ved in four LMC sightlines (not including 
H II regions; iLehner fc Howkii2007i ) . 



able, since any H i lines would be hidden underneath the 
strong H I lines from the neutral phase. In a given com- 
ponent, there are two unknowns, the ionization parameter 
U = n^/nn and the total (neutral plus ionized) hydrogen 
column density A^(H). We ran a grid of CLOUDY models 
at different values of U and Af(H), with \ogU varied from 
—3.4 to 0.0 in 0.2 dex intervals, and logA^(H) varied from 
18.0 to 21.0 in 0.1 dex intervals, and used a two-step iter- 
ative technique to solve for the best-fit values of the two 
variables. In the first step, we found the value of U that 
best reproduces the observed Siiv/Civ ratio or Civ/Ovi 
ratio, depending on which ions are detected. In the second 
step, we found the value of N(R) that reproduces the ob- 
served absolute Civ and SilV (or Civ and Ovi) column 
densities, given the value of U. The best-fit results for CI 
and C2 are shown in Figure 6. 

In CI, the strongest Ovi component where 
iV(Ovi)wl0A^(Civ) and where Siiv is not detected, a PI 
solution can be found at logC/=— 1.8 and logiV(II)=19.5, 
assuming [C/0]=— 0.5. In this model the derived cloud 
density \ognu is w —2.9 and the line-of-sight cloud size 
(A^h/'^h) is «8kpc. This solution, found by matching 
the Civ/Ovi ratio, is consistent with the non-detections 
of Siiv and Nv, and is only possible when the EUV 
enhancement from the background quasar are included. 
UVB-only CLOUDY models do not produce enough O vi, 
as shown by the dotted red line on Figure 6. 

In C2, CiV, Siiv, and Ovi are ah detected. In this 
component the C iv profile more closely follows Si iv than 
O VI. Therefore we use C iv and Si iv as the inputs to the 
CLOUDY model, finding logC/=-2.6±0.1, lognn « -2.1, 



E CompOn 
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<0.046 
0.12±0.07 
0.14±0.04 

<0.33 




Fig. 6. CLOUDY PI models of the high ions in the two 
PDLA components showing Ovi: CI (at — 57kms~^ rel- 
ative to Zabs), and C2 (at — 12kms~"'^). The colored lines 
show the predicted high-ion column densities versus ion- 
ization parameter log U, with the observations (large data 
points) plotted at the best-fit value of logC/. The lower 
X-axis show the corresponding gas density nn- In the CI 
model (upper plot), a solution is found at logt/=— 1.8 
that explains the Ovi, Civ, and upper limits on Siiv and 
Nv. The dotted red line shows the Ovi prediction from a 
UVB-only CLOUDY model (i.e., without the ionizing ra- 
diation from the QSO). In the C2 model (lower plot), no 
single-phase solution is found; a fit to the Civ and Siiv at 
log U—~2.6 under-predicts the O vi column by 0.6 dex (the 
"Ovi shortfall"). These models are described in §8.1. 



and log A^(II)w20.3, but this model under-predicts the ob- 
served O VI column by 0.6 dex (a factor of 4) . This is marked 
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on the lower panel of Figure 6 as the "O vi shortfall" . 
Therefore there is no single-phase PI solution to the high- 
ion columns in this component. Explaining the three ions 
in component C2 by PI alone would require a two-phase 
(or multi-phase) solution, with some fraction of the Civ 
arising with the Siiv, and the remainder arising with the 
O VI. However, such two-phase PI models are ad hoc, non- 
unique, and without a clear physical motivation. Indeed, 
any combination of three ions (in this case O vi, C IV, and 
Si iv) can be reproduced with a two-phase or multi-phase 
PI model, so although we cannot rule out such models, we 
have no clear reason to favor them. 

We do not discuss PI models to C3, C4, and C5 since 
O VI measurements are unavailable in these components, 
and hence the ionization state is poorly constrained. If 
these components are photoionized, the implied gas den- 
sities lognn derived from the Civ/Siiv ratios are —1.9, 
—2.1, and < —2.3, respectively, assuming [C/Si]=— 0.5. 

To summarize this sub-section, PI from the background 
quasar is a viable explanation for the high ions in CI, 
the PDLA component with the strongest Ovi, so long as 
[C/0]=— 0.5 and we include the EUV radiation from the 
background QSO. This model gives a gas density logriH « 
—2.9. No single-phase PI solution is found for C2, since the 
observed O vi is four times stronger in that component than 
the PI model predicts, but the C iv and Si iv in C2 (and 
in C3 and C4) can be explained by PI in plasma at gas 
densities lognn^— 2 if [C/Si]=— 0.5. To explore other pos- 
sibilities, we now turn to coUisional ionization (CI) models. 

8.2. CoUisional ionization models 

In this section several CI models are evaluated in 
their ability to reproduce the high-ion observations in 
the PDLA: single-phase coUisional i onization equilibrium 
(CIE ) and non-e quilibrium mod els (jSutherland fc Dopital 
119931 : iGnat &: S tcrnbcre] 120071), conductive interfaces 
(INT; iBp hringer & Hartquist 1987; Borkowski et al. '19901; 
Gnat et al. 20101'^. and turbulent mixing layers (TML; 
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2006t iKwak fc Sheltonll2010l: iKwak et al.ll201lD . The com- 
mon element uniting these models is that they place the 
high ions in warm-hot, collisionally-ionized plasma rather 
than cool photoionized plasma. In the INT and TML mod- 
els, the high ions trace the radiatively cooling boundary lay- 
ers between cool (^^10^ K) clouds and hot (~10® K) plasma; 
they differ in whether conduction or turbulent mixing is the 
dominant energy transport mechanism from the hot to the 
cool phase. To test the models, we compare the high-ion col- 
umn density ratios they predict with the observations. This 
is shown in Figure 7, with Civ/Ovi vs Nv/Ovi plotted 
on the left and Civ/Ovi vs C iv/Siiv on the right . Thes e 
ratio-ratio plots, first used bv ISavage fc SembachI ()1994[) . 
allow observed ratios to be compared with those measured 
in other environments. 

Before comparing to the observations, the model pre- 
dictions for the high-ion ratios have to be corrected for 
the sub-solar C/0, C/Si, and N/0 abundance ratios typ- 
ical for DLAs, because the models are computed for gas 
with solar relative abundances. When making this correc- 
tion, we neglect any changes in radiative cooling that re- 



^ We use the abbreviation INT to refer to conductive inter- 
faces, since CI is already used to refer to coUisional ionization. 



suit from changing the relative elemental abundances, so 
we assume that the ionic ratios scale linearly with the rel- 
ative abundance of the two element s, e.g. we assum e that 
7V(Civ)/7V(Ovi)oc(C/0), following iFtecet al.l (|2004D . Our 
adopted DLA abundances [C/O]=-0.5, [C/Si]=-0.5, and 
[N/0]=— 1.0 are typical for the neutra l phase of DLAs 
(jPettini et al.l[2008t iPetitiean et al.112008 ^ and are assumed 
to apply in the high-ion phase(s) as well. The effect of these 
abundance corrections on the predicted ionic ratios is shown 
by the solid arrows on Figure 7 (hereafter we refer to these 
new regions of ratio-ratio space as the "corrected predic- 
tions" of the models). 

Our chosen TML models are those presented by 
iKwak fc SheltonI (|2010[ ). We consider the full range of ionic 
ratios predicted by their non-equilibrium Model A (the ref- 
erence model, described in their Table 1). The ratios were 
calculated as a function of time in 1 Myr intervals up to a 
maximum mixing layer age of 80 Myr, and the range shown 
on Figure 7 reflects the variation in the ratios as the mixing 
lay ers evolve w i th tim e. Our chosen INT models are those 
of Gnat et al.l ()2010[ ). We consider the full range of col- 
umn density ratios from both the "Maximal" and "Central" 
cases, covering all parameter choices given in their Table 
4. We show the model predictions as continuous regions of 
ratio-ratio space, rather than discrete points, to account for 
the possibility that the actual values of the model parame- 
ters (coronal gas temperature, density, etc) lie between the 
round-number values chosen when the models were com- 
piled. 

Using the CIE models of iGnat fc Sternberg! (|2003), we 
find there is no single-temperature CIE solution to the 
O VI, C IV, and Si iv column densities, nor is there a single- 
temperature non-equilibrium CI solution in either the iso- 
baric or isochoric cases (not shown on Figure 7). The CIE 
and non-equilibrium CI predictions in ratio-ratio space do 
not overlap with the observed ratios - they miss both the 
PDLA average and the individual component ratios. 

On Figure 7a, inspection of the y-axis shows that the av- 
erage PDLA C iv/0 VI ratio lies within the corrected range 
predicted by the TML models, and toward the upper end 
of the corrected range predicted by the INT models. On the 
x-axis, the PDLA average Nv/Ovi ratio is an upper limit 
since Nv is not detected, but this is consistent with the 
corrected predictions of the TML and INT models. Only 
two individual PDLA components (CI and C2) are shown 
on this figure, since both O VI and N v are undetected in 
C3, C4, and C5 so the ratio Nv/Ovi is undefined in those 
components. 

On Figure 7b, the average PDLA Civ/Siiv ratio also 
lies within the corrected range predicted by the INT and 
TML models, although the ratios in the five individual 
PDLA components show a large scatter around the aver- 
age value. Only two of the five individual components show 
C iv/0 VI and C iv/Siiv predictions that are both matched 
by the corrected INT prediction. The corrected TML model 
fares better, matching both ratios in four of the five PDLA 
components (all apart from CI, which may be photoionized; 
see §8.1). 

The TML model of iKwak fc SheltonI (|2010l) also ex- 
plains the absolute high-ion column densities, at least 
for three of the components: the maximum C iv column 
produced per mixing layer corrected to [C/0]=— 0.5 is 
log A^(Civ)=13.27, and we observe values in the five PDLA 
components of 12.40, 12.69, 12.78, 13.48, and 13.53. The fi- 
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Fig. 7. Diagnostic ratio-ratio plots showing Civ/Ovi vs Nv/Ovi (Fig 7a; left) and Civ/Ovi vs Civ/Siiv (Fig. 7b; 
right), comparing observations (of this PDLA and other galaxy- halo environments) with ionization models. The obser- 
vations shown are the average measurements in the Q0841+129 PDLA (black circles), the component-by-component 
PDLA measuremen ts (open circles), the DLA average from F ox et al. (2007 a ) in b lue, the Galactic disk average from 
iLehner et al.l (|201ll ) in green, the Galactic h a lo av erage from IWakker et al.l ()201l[ ) in orange, and the Galactic high- 
velocity cloud (HVC) average from iFox et al.l (|2005[ ) in red. In each case the error bar shows the st andard deviation o f 
the measurements. T he ionization mod els shown ( shaded regions) are turb ulent mixing layers (Kw ak fc Sheltonll2010[ ). 
conductive interfaces (jGnat et al.ll2010[ ). and CIE fenat fc Sternberg |2007[ ). The shaded regions show the approximate 
range of predicted ratios when allowing the model parameters to vary (see text). The black arrows show the effect of 
correcting the model predictions from a solar relative abundance pattern to a DLA relative abundance pattern with 
[C/O]=-0.5, [C/Si]=-0.5, and N/O]=-1.0. 



nal two columns are marginally higher than the TML model 
can explain. The absolute high- ion columns are h arder to 
explain in the INT models of iGnat et al.l (J201Cl[ ). which 
predict, when corrected to [G/0]=— 0.5, a maximum Civ 
column per interface of log A^(Civ)=12.48, lower than the 
Civ columns observed in all PDLA components except CI. 
Therefore, for the INT model to apply, unresolved multiple 
interfaces would have to exist within each component, in 
order to build up the observed CiV columns. Altogether, 
we conclude that the TML model is more successful. 



It is instructive to compare the model-predicted ionic 
ratios to those meas ured in oth e r DLA s. The average DLA 
high- ion ratios in the lFox et al.l (|2007al ) sample, which con- 
sists of 12 DLAs with Ovi, are shown as large blue points 
on Figure 7. These averages are successfully reproduced 
by both the corrected TML and the corrected INT predic- 
tions, although only just: the TML model cleanly repro- 
duces the average DLA Civ/Ovi ratio (-0.25 in the log), 
but the average DLA C iv/SiiV ratio (0.73 in the log) is at 
the upper end of the TML prediction; conversely, the INT 
model explains the average C iv/Siiv ratio but the average 
C iv/0 VI ratio is at the upper end of the INT prediction. 
Given the uncertainties involved in correcting the model 
predictions to the non-solar DLA abundance pattern, we 
judge these models to be fairly successful at reproducing 
the average DLA ionic ratios. 



Based on the above discussion, we conclude that bound- 
ary layers are a viable explanation for the C iv and O vi 
absorption, both in this PDLA and in DLAs in general, 
and that turbulent mixing is the most successful of the 
models examined here, explaining the high-ion ratios in 
four of the five PDLA components. The boundary-layer 
explanation is supported by the finding that the average 
Civ/Ovi and Civ/Siiv ratios observed in DLAs are sim- 
ilar to those measured in the disk, halo, and high-velocity 
clouds (HVCs) of the Milky Way (a DLA itself), where in- 
terfaces and turbul ent mixing are widely invo ked to explain 
the high ions (e.g. ISayage et al.|[l9 94. 2003 llSpitzeiJ 119961: 
'Sembach et al.' '2003; 'Colhns et al. 2004, 2005'; 'Fox^.e£aD 
2004, 2005; Lchncr et al. 2011; Wakkcr et al. 2011]), al- 
though inspection of Figure 7b shows that the Galactic 
disk and halo average C iv/Si iv ratios are marginally lower 
than the predictions of the solar-abundance TML and INT 
models shown. The DLA average Nv/Ovi ratio is at least 
0.5 dex lower than the Galactic disk and halo average 
Nv/Ovi ratios, but this can be explained by the sub-solar 
N/0 ratio prese nt in DLAs, i.e. as an abundance effect 
(|Fox et al.ll2009h . 

The boundary-layer scenario, in which the high ions 
are cooling rather than photoionized, has other advan- 
tages: it explains the metallicity-indep endence of the 
stren g th of O vi absorption in galaxy halos (jHeckman et al.l 
120021: Wo^ 1201 ID . and fits naturally into the established 
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picture of galaxy forniation where cold streams pass 
through hot plasma, since boundary layers will naturally 
arise between the two phases. Indeed, in this picture, 
the Ovi and CiV detections in DLAs indicate (without 
tracing directly) the presence of a hot (~10^K) phase 
of surrounding plasma in the host galaxies. Such hot 
galactic coronae are predicted to arise by shock heat- 
ing of gas falling into poten t ial wells (hot-mode ac- 
cretion: iRees fc Ostrike^ Il977t iBirnboim fc Dekell l2003t 
iKeres et all I2005L 12001 ) . and by feedback from stel- 
lar activity (supernoyae) within the galaxies themselves 



e.g. Rasmussen et al.ll2009l : ICrain et al.ll201Cll : lHenleY et al.l 



l2010f ). If this boundary layer scenario is correct, then the 
larger significance of Ovi observations in DLAs at z=2- 
3 would then be that feedback and/or hot-mode accretion 
are already active at this epoch and have created hot galac- 
tic coronae. Further studies of intervening DLAs where the 
high- ion ratios can be measured component- by-component, 
and where the ionizing radiation field and the influence 
of photoionization are expected to be lower, will confirm 
whether this picture is correct. 

Finally, we note that neither the boundary-layer sce- 
nario nor the photoionization scenario explains the high-ion 
kinematics, which are considerably more extended than the 
low-i on kinematic s, both in the PDLA and in DL As in gen- 
eral ()Wolfe fc Prochaska 2000al: iFox et al.|[2007b[ ). This in- 
dicates that an energetic mechanism (e.g. a galactic wind or 
galaxy-galaxy encounter) has dispersed the high ions over a 
large velocity interval. Simulations of the high-ion velocity 
fields in DLAs are needed to investigate the nature of this 
mechanism. 



9. Summary 

In order to explore the processes that generate high-ion 
plasma in high-redshift galaxy halos, we have presented 
new VLT/UVES high-resolution observations of the PDLA 
at 2;abs=2. 47621 toward Q084H-129, focusing on the origin 
of the high-ion absorption. This system lies at 1 630 kms^"'^ 
from the background QSO, has logiV(Hi)=20.80±0.10, a 
metallicity [Zn/II]=— 1.60±0.10, a low-ion velocity width 
A?;go=30 kms^^, and a hydrogen ionization fraction in the 
low-ion phase of 0.1-0.5 as determined from the An/Si ii, 
Ar i/S II, N ii/N I and Fe iii/Fe ii line ratios. The high ions in 
the PDLA show a multi-component structure; at least five 
C IV components are seen spanning a total velocity extent 
of ssl60kms~^ with 6-values ranging from 8.8±1.8kms~^ 
to 27.3±3.1kms~"'^. Siiv is detected in three of these com- 
ponents; Ovi is detected in two, and blended in the others. 
The PDLA is (partly) shielded from the QSO's ionizing 
radiation by an intrinsic sub-DLA at Zabs=2. 50620, 
very close to the QSO emission-line redshift of 
Zom=2.49510±0. 00003. This filtering arrangement hardens 
the ionizing radiation field incident on the PDLA, allowing 
us to test whether photoionization (PI) is a viable origin 
mechanism. We find that a single-phase PI model at 
logriH ~ —2.9 can successfully explain one of the two com- 
ponents showing Ovi (CI), so long as [C/0]=— 0.5, but 
no single-phase PI model is found for the other component 
showing Ovi (C2), since the Ovi in this component is a 
factor of four stronger than the model predicts. The Civ 
and Si IV profiles in C2, C3, and C4, are consistent with 
PI at logriH ~ —2 if [C/Si]=— 0.5. Single-phase collisional 



ionization models (equilibrium or non-equilibrium) are 
unable to reproduce the high-ion observations. 

An alternative model, where the high ions trace the tur- 
bulent mixing lay ers between the warm D LA gas and a hot 
exterior medium (jKwak &: Sheltonll2010t) . can successfully 
reproduce the ionic ratios in four of the five components 
(C2, C3, C4, and C5). This model can also explain the av- 
erage high-ion ratios Civ/Ovi, Nv/Ovi, and Civ/Siiv 
observed in a larger sample of 12 DLAs. Turbulent mixing 
layers are often invoked (along with conductive interfaces) 
to explain high-ion observations in the halo of the Milky 
Way (itself a DLA). If correct, this model provides indirect 
evidence for the existence of a hot phase of coronal plasma 
in and around DLA galaxies. Further observations of high- 
ion absorption in individual intervening DLAs are needed 
to test this model. 
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